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Abstract 

This is the third in a series of papers in which we set up and discuss the functional 
relations for the "split rapidity line" correlation function in the TV-state chiral Potts 
model. The order parameters of the model can be obtained from this function. Here 
we consider the case N = 3 and write the equations explicitly in terms of the hyper- 
elliptic functions parametrization. We also present four-term low-temperature series 
expansions, which we hope will cast light on the analyticity properties needed to solve 
the relations. The problem remains unsolved, but we hope that this will prove to be a 
step in the right direction. 

PACS: 05.50.+q, 02.30.Gp Keywords: Lattice models, chiral Potts model, order pa- 
rameters, hyperelliptic functions 

1 Introduction 

In a previous paper we used the method pioneered by Jimbo, Miwa and Nakayashiki 
]|] to obtain functional relations for a generalized one-spin correlation function for the 
planar solvable iV-state chiral Potts model. The correct solution of this would give the 

1 



spontaneous magnetization - a long-standing problem which is all the more tantalizing 
as there is an elegant conjecture for the result Q - 0. 

Unfortunately, the functional relations are really just symmetry and periocity con- 
ditions and have more than one solution. One needs to supplement them with an 
appropriate analyticity requirement. For other models this is easy: one goes to the 
elliptic (or trigonometric) parametrization that uniformizes the equations and mani- 
fests the rapidity-difference-property. There is then an obvious vertical strip in the 
complex plane of the spectral parameter variable within which the function is ana- 
lytic and Fourier transformable in the vertical direction. The equations can then be 
solved directly and uniquely by taking such a transform, in much the same way as the 
inversion relations for the free energy can be solved Q. 

For the chiral Potts model with N > 2, there is no such parametrization. In one 
sense one can uniformize the equations, but at the price of introducing hyperelliptic 
functions with more than one argument, the arguments being related to one another 
in a complicated manner. We have discussed these functions (particularly for the case 



N = 3) earlier |Q] |1C], and in [11] have obtained formulae for the coefficients in 
the functional relations. Here we write the relations in terms of this parametrization, 
focussing on the case N = 3. We discuss their properties and present some low- 
temperature expansions of the functions, obtained by direct calculation from their 
definitions. We still have not succeeded in obtaining the appropriate solution of the 
functional relations, but hope that the data presented herein will assist in the search. 
Certainly the solution must agree with the low-temperature expansions. 



2 Functional relations 

The chiral Potts model is explained in |lj. Spins live on the sites of a square lattice, 
oriented diagonally. Each spin takes the values 0, . . . , N— 1. Adjacent spins a, b interact 
with weight function W pq (a — b) on SW — * NE edges, W pq (a — b) on SW — > NE edges. 
The symbol p denotes a point (a "rapidity") a p ,bp,c p ,dp on the projective curve 

a^ + k'b^ = kd% , fc'< + ^ = kc£ , (1) 
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where k, k' are real constants (moduli) satisfying 

k 2 + k' 2 = 1 . (2) 
The symbol q denotes another such point and 

«w») = n 't'^ ■ w„ ( n) = w (3) 

Operators R and M act on the rapidities, being defined by 

(aRp,bRp,CRp,dRp) = (b p ,ujap,dp,Cp) , (4) 

(a>Mp,bMp,CMp,dMp) = (ua p ,bp,c p ,ujdp) . (5) 
They satisfy RM = M~ 1 R and ensure that 

W pq (n) = W q , Rp (n) , W pq {n) = W q , Rp (-n) . (6) 

In Q we consider a square lattice with one spin deep inside it fixed at some value 
a. The horizontal rapidity line immediately beneath this spin is broken, the left and 
right segments having different rapidities p, q, respectively. The boundary spins are 
fixed to zero. Writing the corresponding partition function as Z pq (a), the probability 
that a free spin in this position has value a is simply 

F pq (a) = Z pq (a)/[Zp q (0) + ■■■ + Z pq (N - 1)] . (7) 

In the limit of a large lattice, it is independent of the rapidities of all other lines 
(because of Z-invariance By definition, F pq (0) + • • • + F pq {N — 1) = 1. 

In M we show that the Yang-Baxter relations imply that F pq (a) satisfies the fol- 
lowing functional relations: 

FR P ,Rq(a) = Fpq(-a) , F Rq ^(a) = F RPjq (a) , 

N-l 

F P M a ) = ^E 1 ^-^^) - (8) 

6=0 

F Mp ,q(a) = F pM -i q (a) = r] pq co a F pq (a) , 



where 



N-l N-l 
i pq = zZ W Pi( n ) > Vpg = Y.^FpM • ( 9 ) 

n=0 a=0 



Let u! = exp(2iri/N). It is natural to introduce the discrete Fourier transforms (for 
r any integer) 

N-i 

F pq (r) = Y." ra F pq {a) , 

a=0 
N-l 

Vpiir) = $> ra W M (a) ■ ( 10 ) 

a=0 

Then F pg (0) = F pq (N) = 1 and the relations (|8|) become 

F Rp , Rq {r) = F pq {N-r) , F^^r) = F Rp>q (r) , 

F p , Rq (r) = V pq (r)F q , Rp (r)/V pq (0) , (11) 
F Mp , q {r) = F PjM -i q (r) = F pq (r + 1)/F pq (l) . 

Without loss of generality, in these equations we can take r = 1, . . . , N — 1. 
It can be convenient to work with the ratios: 

G Pq (r) = F pq (r)/F pq (r-l) (12) 

and with the product 

L pq{ r ) = G pq {r)G RqtRp {r) . (13) 

In terms of L pq (r), the functional relations become, for r = 1, . . . , N, with 
L pq {N + 1) = L pq (l): 

L Rp , Rq {r) = L qp (r) = l/L pq (N-r + l) , 

L p , Rq {r) = V pq {r) L q:Rp (r)/V pq (r - 1) , 

L pq (l) ■■■ L pq (N) = 1 , (14) 

LM P , q ( r ) = L P ,M-igW = L pq (r + 1) . 



It follows that 



L R^P, q {f) = L P)R ?N q (r) = (-1) N 1 L pq (r) , (15) 



which is a periodicity (or anti-periodicity) property. 

The function F pp (a) is the probability, for the regular square lattice, that a central 
spin has value a. This should be independent of the rapidity p fl^l , so F pp (a) , F pp (r) , 



4 



Gpp(r), Lpp(r) should all be independent of p. The order parameters are the ex- 
pectation values M r of u ra for r = 1,...,N — 1. From the above definitions of 
F pq (r),G pq (r),L pq (r), we see that 

Mr = F pp (r) = Gpp(l) • • • G pp (r) = ^ L PP {1) ■ ■ ■ L pp (r) . (16) 

The conjectured result for M r is 

Mr = <u/ a > = (1 - k'Y N - r ^ 2N2 , (17) 

(eqn 3.13 of || , eqn. 1.20 of [f| , eqn. 15 of ||, f3 and A therein being the k' of this 
paper. This k! can be thought of as a temperature variable: small at low temperatures 
and increasing to one at criticality. 



3 Hyperelliptic function parametrization 

The hyperelliptic function parametrization is given in @ and 0. There are N — 1 



constants p±, ... , pn-i (positive pure imaginary) which are defined solely by N, k and 
k', and satisfy p a = pN~a- For a, f3 = 1, . . . , N — 1, define 

Ta(3 = Pa + P/3 ~ P\a-/3\ , (18) 

taking po = Pn = 0. Let s = {s±, . . . , sat_i} be a set of variables. Then the hyperel- 
liptic function is defined by 

= ^2exp{2TTi^2m a s a + TTi^2^2m a T al3 mp} , (19) 

m a a R 

the inner sums being over a, = 1, . . . ,N — 1, and the outer sum over all values of the 
integers m = {mi, . . . , niN-i}- 

We also define the constant sets g,p by g = {g%, . . . ,gN~i}, P = {pi, ■ ■ ■ iPn-i}, 
where g a = a/N. From (18) and ([l9|) the @ function satisfies various quasi-periodicity 



and symmetry relations, in particular 
9(s) = e(-a) = 0(5) = e{s + Ng) = exp{mJ2(. 2s a + Np a )} @(s + Np) , (20) 

a 

where s = {sjv-i, sjv-2j • • • ; si} is the sets in reverse order. 



For each rapidity p there is a set of such s-variables, which we write as s p (Greek, 
numerical or "iV" suffixes refer to the position of an entry in the set, while lower-case 
Roman suffixes - usually p or q - refer to the rapidity dependence). We write the 
entry a of s p as (s p ) a . The entries (s p )i, • • • , (s p )n-i are not independent: they are 
functions of the single rapidity variable p and satisfy the N — 2 relations (36) of P], or 



equivalently (10) of 



Set t = {ti, . . . , tjv_i}, where 

t a = ai + • • • + s a - a(si H h a N -i)/N . (21) 

As with s, we write t p for the set t corresponding to the rapidity p, and {t p ) a for its 
entry a. Then 

(tRp)a = (tp)jV-a + \Pa , (22) 

^r 2 p = tp + P i ^Mp = t p - g . (23) 
(choosing M so that M~ l is the operator M x (1) • • • of and |ll[.) 

The relation (|22|) , together with the s, s symmetry in (|20|), implies that 

©(^g - + mfif + n/)) = @(t q - t p - mg + np) (24) 

for all integers m, n. 



From pi) and (TO) 



Vpq{r)/V pq (r - 1) = (C p 6q - LO r dpdq) / {bpC q - uf d p d q ) 

so 

Vp,Rq(r)/V p , Rq (r - 1) = (uCpdg - uj t a p Cq) / (b p dq - iv r d p bq) . (25) 



In eqn (26) of |11] we give a formula for the RHS of (25). This formula has not been 
proven for N > 3, but we have tested it numerically and believe it should be proveable 
by Liouville's theorem, along the lines indicated in Using it, we obtain 

Vp,Rq{r)/V P)Rq {r - 1) = ip P ,R2q(r)/tp pq (r) , (26) 

where 



VvA) = nr , Kq , .V ^ 1 • (27) 
0[* 9 -t p + (r- 1)5] 



0(tg - t p + rg) 

[tq ~ t p + (r - 1 

Replacing g by the second line of equation (|i~4|) is (using the first equation) 

L P ,R2q{r) = ^p,R2q(r)L pq (r)/ijpq(r) . (28) 



A spurious solution 

Obviously a solution of fl28| ) is 

L Pq (r) = iP pq (r) . (29) 

Further, using (f23|) and (|24|), we find that it satisfies all the other equations (|l4|). For 
N = 2, it is indeed the correct solution, but for general iV the formula (25) of |nj (also 
unproven for N > 3, but believed to be "proveable") tells us that 

^__(r) = fc( JV_1 )/ 7V TT MiZ^ cp°g-^ %£i . (30) 

MV 1 Cpd, - uJ- m ~ l d p c q b p a q - uji- m - l a p b q J 

This is precisely the wrong solution discussed in jjj]. It gives a result for M r that is 
in error by a power iV/2, and disagrees with low-temperature series expansions. Yet 
it is such a simple and elegant (in terms of our hyperelliptic functions) solution of 
(|j~4|). This illustrates well the difficulty with functional relations such as (0): they do 
not by themselves define the function. One has to incorporate the correct analyticity 
properties. 



4 The case N = 3 

If N = 3, then pi = P2 = p, so if we define 

x = e 27Tip , \x\ < 1 , (31) 

then 

0{ Sl ,s 2 } = $(6^,^) , (32) 

where 

= ^2;™ 2 +»™+™ 2 a ™/3™ j (33) 



the sum being over all integers m, n. (Here we regard x as a given constant and do 
not usually explicitly exhibit the dependence of functions on it.) The function &(a,(3) 



satisfies the symmetry and quasi-periodicity properties (34) -(36) of [11], in particular 



$(a,P) = *(l/a,l//3) = *(a,a//3) . (34) 
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The nome x is related to the modulus k by 



1 12 



{k'/kf = Tlx Q{x 6 )/Q{x) , (35) 
where 

oo 

Q(x) = ]J(l-x n ) . (36) 

71=1 

At low temperatures both k and x are small. 
If we define 



(37) 



a = 




= exp{27ri[(sg)i 


~ (s P )i}} , 




a = 


fipq 


= exp{27ri[(sg)2 


" (*p)2]} , 




u = 




= exp{27Ti[(s g )i 


- 2(s q ) 2 - (s p )i + 2(s p 


) 2 ]/3} 


z) = 


oo 

71=1 


-x 3n ~ 2 z)(l -x 3n 


- l /z) , <p(z) = p(z- 





(38) 



then 



a = u 3 [3 2 , (39) 

W pq {\) = n 2 ^0(a)^(a//3) , W pq (2) = uf3 <t>{a)<j>{fi) . (40) 

For any function, we can in principle eliminate the two degrees of freedom p and q 
in favour of the variables a, (3 (regarding k,k',p,x as given constants). The resulting 
expressions will not necessarily be simple: they may be multi-valued functions of a 
and f3, but the simple form of ( |40| ) gives some encouragement that this may be a useful 
parametrization for our functions. In this spirit, let us define a function ip{a,(3) such 
that 

F pq (l)/F pq (0) = xu~ l [3- l ilj{a,P) . (41) 

We have expanded tp{a,(3) to third order in the low-temperature parameter x, 
working directly from the definition (0) and performing finite-size lattice calculations. 
The results are given below. We used a Fortran computer program, working with series 
in x with double-precision numerical coefficients for various given values of a and (3. 
The method is not rigorous, depending both on observing that a given coefficient 
stabilized once the lattice was sufficiently large (the biggest lattice we considered was 
about eight by eight), and on numerically fitting the results to postulated expressions 
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with unknown integer coefficients. However, we believe the results to be correct. The 
key observation is that for n > the coefficient of x n in the expansion is a Laurent 
polynomial in a and j3, divided by (a — l) 2n_1 . The coefficient of x° is one. 
The variables a pg , (3 pq , u pg satisfy the relations 

&Rp,Rq — &pq i @Rp,Rq — ^pq/ Ppq i URp,Rq = l/Upg , 
OLqp — \/(Xpq , Pqp — 1 / Ppq , Uqp — j 
a p,R 2 q = x a pq : Pp,R 2 q = x Ppq ; u p,R 2 q = u pq ; (^2) 
&Mp,q = Olp^M-lq = a pq , Phlp,q = Pp,M~ 1 q = Ppq i 
UMp,q = u p,M~ x q = w_1 u pq ■ 

Using these, the first of the equations @ is equivalent to 

Fpq(2)/F pq (0) = xu- 2 (3- l il){(x,a/P) . (43) 



Using fl42|) ) (4C) and (||), the second and third of the relations (JSj) become 

xaip(x~ 2 /a, x~ 1 (3/a) = ip(a,(3) , (44) 

A(a,P)ip(x 2 /a,x/(3) = </)(a)(j)((3) + (x/a) tp{a, a//3) + (x//3)(/)(a)(j)(a//3)ip(a, P) , 

(45) 

where 

A(a,p) = l + xif>(a)<l>(a/P)$(a,a/p)+x<i>(a)<l>(J3)il}(a,P) . (46) 

The last two of the relations (||) (those involving Mp and M~ 1 q) are automatically 
satisfied. 

Series expansions for ijj(a,(3) 

Define 

, n = pvv . (47) 



P{a - 1) /?(<* - 1) 
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Then, using the configuration of the rapidity lines p and q mentioned above, taking a 
and P to be of order one, to third order in an expansion in powers of x we find that: 

t(}(a, f3) = 1 + x[va(/3 - 1) - 2 + a] + x 2 [u/x(a 2 - 0) + v(l + 4a + 2a 2 - 
2(3 - ba/3) + 13 - 6a] + x 3 [v^ 2 a(3(3 - a/3 - a - a 2 ) + f/x(2/3/a + 5/? + 
3a/3 + 1 - 2a - 8a 2 - a 3 ) + u(-6/9/a + 12/? + 28a/3 + a 2 /3 - 3 - 21a - 
10a 2 - a 3 ) + 5/a - 71 + 38a + a 2 ] + 0(x 4 ) . (48) 

We also considered the configurations where one of the lines p, q was rotated through 
180°, as discussed in jL|. This gives results for four cases where the arguments of the 
function tp are themselves proportional to non-zero powers of x, namely: 



<ip(xa,x(3) = l + x[(l-a)//3-2]+x 2 [(l-a)//3 2 + (-4 + 6a- 
a 2 )//3 + 1/a + 9 + 3a] + x 3 [(1 - a)//? 3 + (-5 + 7a - a 2 - a 3 )//? 2 + 
(2/a + 17 - 32a + 11a 2 - a 3 )//3 - 10/a - 49 - 23a + 2a 2 + (2/a + 4 - 
a)/?]+0(x 4 ) , (49) 

tp(xa,P) = l-2x + x 2 [(l-a)//3 + l/a + 9 + 3a + (2/a + 4-a)/?- 
(1 - a) 2 /3 2 /a 2 ] + x 3 [(-4 + 6a - a 2 )/ (3 - 10/a - 49 - 23a + 2a 2 + 
(3/a 2 - 19/a - 29 + 14a - a 2 )(3 + (-1/a 3 + 13/a 2 - 20/a + 11 - 
3a)/3 2 - 2(1 - a) 2 /3 3 /a 3 ] + 0(x 4 ) , (50) 

ip(a/x,/3) = 1 + a + x [(1 - a) 2 / [3 - 4 - 4a + 0\ + 

x 2 [(1 - a)(l - a 2 )//3 2 + (1/a - 8 + 13a - 8a 2 + a 3 )/f3 - 2/a + 25 + 

25a - 2a 2 + (1/a - 5 + a)0\ + x 3 [(1 - a)(l - a 3 )//? 3 + (3/a - 14 + 

11a + 11a 2 - 14a 3 + 3a 4 )/^ 2 + (1/a 2 - 14/a + 58 - 79a + 58a 2 - 

14a 3 + a 4 )//3 - 2/a 2 + 19/a - 146 - 146a + 19a 2 - 2a 3 + 

(1/a 2 - 8/a + 30 - 8a + a 2 )/3] + 0(x 4 ) , (51) 

Tp(a/x,/3/x) = 1 + a + (3 + x [-4 - 4a + (1/a - 5 + a)(3] + 
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x 2 [(1 - a) 2 / (3 - 2/q + 25 + 25a - 2a 2 + (1/a 2 - 8/a + 30 - 8a + 

a 2 )/3] + x 3 [(1/a - 8 + 13a - 8a 2 + a 3 )//? - 2/a 2 + 19/a - 146 - 

146a + 19a 2 - 2a 3 + (1/a 3 - 11/a 2 + 56/a - 175 + 56a - 11a 2 + 

a 3 )/? + (1/a 3 - 2/a 2 + 2/a + 2 - 2a + a 2 )/? 2 + (-1/a 3 + 3/a 2 - 

4/a + 3 - a)(3 3 } + 0(x 4 ) . (52) 



In all of the expansions ( |48| ) - (pig), the variables a, (3 are of order one. 

We can regard these expressions as giving the expansion of the function vp(a, (3) in 
various domains in the (a, (3) plane. We have verified that the results for neighbouring 
domains are consistent with one another. For example ( |49| ) and (^) are consistent in 
that if one replaces (3 in ( |50j ) by x(3, then one obtains (p9|), except for terms in ( |49| ) 
that are of order higher than 3 when (3 is of order a? -1 . Similarly, (fig ) is consistent 
with both (|0|) and (ID, and (||) is consistent with (||). 

Taking a = 0(x _1 ) and /3 = 0(1) or 0(x _1 ) in ((44|), one can readily check that 
( ]5l| ) and ( |5"2"D each has the consequent symmetry property. Also, taking a = O(x) and 
(3 = 0(1) or 0(x) in (|^), we find that the relation is satisfied by (p9|) and (|50|). 

If q — > p, the p variables remaining finite, then a, (3 — > 1 while u and U remain finite. 
The coefficients of the terms proportional to v and U in ( |48| ) then vanish, so we obtain 
the unique result: 

^(1,1) = 1 - x + 7x 2 - 27x 3 + 0(x 4 ) . (53) 

From (fn|) and (|43|) it follows that F pp (l)/Fpp(0) = F pp (2) / F pp (0) = xi/)(l, 1), so from 
(H) 

^1 = M * = TT^mi^ = l- 3x + 9x 2 - 45x 3 + 231x 4 + 0(x 5 ) , (54) 



in agreement |8[ with the conjecture (|17|) . 

Fourier transformed relations 

To obtain the relations (|TT|) , ( |i~4j) in terms of the variables a,/3 - or rather in terms of 
u, (3 - define: 

A(u,P) = l + xf3- l u- l tp(a,(3)+x(3~ 1 u- 2 ij(a,a/f3) , (55) 
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B(u,/3) = A(u,P)A(u,u- 3 p- 1 ) , 
({u,(3) = [l + uf3 (j){a) 4>(fi) + u 2 (3(/)(a)0(a/f3)} /A(a,/3) 
ThenF pg (0) = l/A(u,P), 

F pq (r) = A(u;- r u,P)/A(u,P) , 



L pq (r) = B(tu~ r u, (3)/B(u u, (3) 



and, from and (I 



A(l/u,u 3 /3) = A{l/u,x~ l (3- 1 ) = A{u,p) , 
A(l/u,x/f3) = ((u,P)A(u,p) , 



(56) 
(57) 

(58) 
(59) 



(60) 



B{u' x ,u 3 p) = B^u" 3 /?" 1 ) = B(u,(3) , 
B^^xp' 1 ) = ({u,P)B( U - 1 ,p- 1 ) . 



(61) 



Together with (42) and the fact that A(a, P) is a single valued function of u 3 and P, 
the relations (^) imply (pi]), while ( |6l"D imply (|14|). One does not need the definition 
of A(a,/3). 



Our function p{z) is the ip(z) of |11]. Another function discussed therein is 



h(a,P) = ${u,u z P)<&{uju,uj 2 u 2 P)§{uj' i u,uju z P) . 



(62) 



Using also eqns (32) and (58) of [11], it follows that 



C(u,P) = T(a,P)$(u,u 2 p/x)/<S>(u,u 2 p) 



where 



r(a,P) 



h(a, P) 



Q{x) 2 Q(x 3 ) 4 p{a) p{p) p(a/P) A(a, P) 



(63) 



(64) 



Obviously (|58|) and (|59|) are unaffected by multiplying the functions A(u, P), B(u, P) 
by factors which are single-valued functions of « 3 ,/3 (i.e. of a,P ), such as h(a, P) and 
r(a, P). To within such factors, using the symmetry property (|34[), we see that a simple 
solution of (|l|) is B(u,P) = ®(u,u 2 P). This implies that L pq (r) = $(uj~ r u,uj r u 2 P)/ 
^(io 1 ~ r u,uj r ~ 1 u 2 P), which is the spurious solution ([E]) mentioned above. 
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More explicitly, let us define 

C(u,p) = — 5—- , D(u,p) = T-pr . (65) 

$(u,u 2 /3) ,y ' [C{u,P)C(uju,P)C{lj 2 u,P)} 1/3 

Thus D(u,P) is C(u,/3) normalized by the multiplication of a function of a,P so that 

D(u, P)D(tou, P)D(lj 2 u, ft) = 1. The equations ( |61~| ) then reduce precisely (including 

all factors) to the simple form 

D(l/u,u 3 p) = D{u : u- 3 p~ l ) = D(u,P) = D{u,xP) . (66) 

From our series expansion (|48|) , for u,a,P of order 1, we find to order x 4 that 

D(u, 1) = 1 - {u + u -1 )x + (l+2n+2u" 1 )x 2 + (u 3 /3 + ^/S-lh-llu- 1 -4) x 3 + 
(-4u 4 /3 - 4u~ 4 /3 - 2u 3 - 2u~ 3 - Wu 2 /3 - 10u~ 2 /3 + 56u + 56u _1 + 26) x 4 (67) 

and 

D(u,p)/D(u, 1) = 1 + [(a+l-2p)u + (a+l-2a/ p)/u] /xx 3 + [(2a/3 + 2/3 - 4a)u/j 2 
+(2a 2 //3 + 2a//3 - ia)fi 2 /u + (14/3 - 8a - 8 + p 2 + /3 2 /a)u/u + (14a//3 - 8a - 8 + 
a 2 IP 2 + a/p 2 )n/u + (1 - /3)(a - /3)(u//3 2 + 3u/a + a" 1 /^ + 3/3~ 2 /u) ] x 4 , (68) 

while from (49) - (|52"D, for it, a,/3 of orders x 1/,3 ,x, 1, respectively, we find to order x 4 



that 



D(u, 1) = l-(ii + u" 1 )x + (l+2n+2n" 1 )x 2 + (n 3 /3 + n" 3 /3-ll?i-ll^ 1 -4)x 3 + 
(-4u~ 4 /3 - 2tT 3 - 10u" 2 /3 + 56u _1 + 26) x 4 + (2u~ 5 + 47u~ 4 /3 + 58u~ 3 /3)x 5 - 
{3u- 7 + 7u- 6 /9) x 6 



and 



D(u,P)/D(u, 1) = 1 + (1 - /3)(1 - /3x/a) {(«- 2/3u + u" 1 - 3cm -1 //? + 3au" 1 )x 3 + 
(-P 2 u/a-2pu/a-Wu- 1 +3PU- 1 -3u- 1 /p)x i + 3{l + p + p 2 )u~ 1 x 5 /a} . (70) 



Applying the symmetries ( |66|) to these expansions, we can obtain D(u,P) to order x 4 
for u of order 1, x 1 / 3 and x -1 / 3 , and for all /3. We have verified that the expansions 
are consistent with these symmetries, and the results for neighbouring domains are 
consistent with one another in the sense discussed above, provided we remember that 
some coefficients have poles when a = 1, x 2 , x -2 , .... 
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Order parameters in terms of D(u,f3) 

From (16) and the above equations we find that the order parameters Mi,M.2 are 
given in terms of the function D(u,f3) by the simple relation 

M\ = Ml = , (71) 

the factor A; 2 / 3 being the contribution from the spurious result obtained by taking 



D(u,(3) = 1, i.e. by using (29). The right-hand side is to be evaluated in the limit 
q — > p, i.e. when the arguments u,/3 of D(u, j3) behave so that u 3 ,/3 — > 1, the ratio 
(/? — l)/(ii 3 /3 2 — 1) remaining finite. 

Using the expansion (|67[), we again obtain (|54[). 



5 Discussion 

If we can solve the functional relations for the "split rapidity line" correlation function 
F pq (a), then we can calculate the order parameters Mj- For the N = 2 chiral Potts 
model, i.e. the Ising model, this can be done straightforwardly, using a uniformizing 
elliptic function parametrization. For N > 3 it is still an unsolved problem. For 
N = 3 we have written down the equations explicitly using a hyperelliptic function 
parametrization: in (^4|) - (pq) in terms of the function ip(a,(3), and in (pq ) in terms 
of D(u,P). We have also presented the first few terms in a series expansion, obtained 
from finite lattice calculations. 

The problem is how to solve the equations. One modest approach is to try to 
guess the general form of the coefficients in the series expansions, and then see if 
the relations determine the coefficients precisely, at least successively term-by-term. 
Such an approach works well for the Ising model in terms of its original variables (the 
Boltzmann weights) ]14| , eqns 2.15 to 2.30]. As a first step, we tried replacing every 
numerical coefficient in (|4g| ) - (^) by an arbitrary number, substituted the expansions 
into (|i||) - (|46"1), and attempted to systematically solve the resulting equations for the 
unknown numbers. This works for ip(a,/3) to zero and first order in x, but at second 
order we were left with four undetermined coefficients, making it impossible to proceed 
further. (We did try looking ahead to see if the higher-order equations fixed these four 
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unknowns, but with no success.) 

This approach appears to be even less successful when applied to the relations ( 66 ): 
for D(u,/3) to first order, which corresponds to zero order for i/)(a, (3), we easily see 
that a solution is D(u,(3) = 1 — fj,(u + u~ 1 )x, where \x is an arbitrary constant. In 
fact it is painfully obvious that the equations are unaffected by multiplying D(u, (3) by 
any function f(u) satisfying f(u) = f(u~ 1 ), f(u) f(uju) f{uj 2 u) = 1. This can certainly 
affect M\ and M.2- 

It could be very useful to merely know the functions to leading order for all values 
of a, j3, or u, (3: this would provide valuable clues as to the locations of any zeros 
or poles. This would be very significant if the functions are single- valued functions 
of their arguments, but we have no reason to suppose this, and know of no way of 
obtaining even the leading-order behaviour except near the "physical regime" where 
the Boltzmann weights are real and positive (and those regimes which map to it by 
the symmetries used to derive the functional relations). 

It may not be necessary to consider the full complex u and (3 planes. If we discard 
the first equality in (^), we have two equations in which we can regard u as a fixed 
"constant". They relate the values of D with second argument j3, x(3 and ii _3 //3, 
so are much like inversion relations ||, the inversion points being f3 = u~ 3 / 2 and 
(3 = x 1//2 u -3 / 2 . If D were analytic in an annulus in the complex (3 plane, centre the 
origin, including these points, then it would be easy to solve (|66|) using a Laurent 
expansion in this annulus. In fact the solution would be a constant, i.e. D(u,(3) would 
be independent of (3. 

To second order in x this is the case, but from ( |68| ) we see that at third order 
there are terms proportional to ^, which from ( f47| ) and ( |3^ ) has a second-order pole at 
(3 = u -3 / 2 . (At fourth order there are terms proportional to /j?.) So D is not analytic 
at the inversion points and does depend on (3. 

We are left with a tantalizing puzzle: undoubtedly the functional relations do con- 
tain information, but they need to be supplemented with a knowledge of the analyticity 
properties at the "inversion point" . Perhaps we should be using the individual hyper- 
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elliptic function variables |J 

z p = x 1//2 exp[27ri(s p )i] , w p = exp[27ri(s p )2] , (72) 

in terms of which a = z q /z p , (3 = w q /w p . They can be chosen to be of order unity 
when x is small, and to leading order w p = z p + 1 = (a — I) /(a — /?), w q = z q + 1 = 
j3(a — l)/(ci! — 0). Higher-order terms can be obtained from equations (4.5) and (4.6) 
of I- 

A third variable that enters the hyperelliptic parametrization of the Boltzmann 
weights is 7 = w p w q /z p [11]. The variables a,P, 7 all have the property that they 



are unchanged by simultaneously replacing p, q by R 2 p, R 2 q, so they automatically 
incorporate this symmetry of the generalized correlation function. One can write the 
expansion in ( |48| ) (at least to the order given) in a form where the coefficients are 
Laurent polynomials in a, P, 7. The trouble is that this expansion is not unique, since 
a,P, 7 are related to one another. Defining \i as in (|47|), i.e. \x = P~ l {P — l)(a — 
P)/(a - l) 2 , we find that 

H = 7" 1 + [7/a - 3 - 3a" 1 + (2a + 2a" 1 - l)/7 + (1 + a)/7 2 ] x + 0(x 2 ) . 

Of course it may be that this hyperelliptic parametrization is not helpful at all, but 
this seems unduly pessimistic since (a) the parametrization does provide convenient 
and simple expressions for the relevant Botzmann weight functions and their Fourier 
transforms [0], and (b) the coefficients in the expansions are much simpler than they 
are if one uses the original variables a p , . . . , d p , a q , . . . , d q . While it is true that the 
coefficients in our series expansions contain negative powers of a — 1, there do not 
appear to be any negative powers of P — 1 or a — P (which occur in, say, the expansion 
of 7). It also appears from ([37]) that such negative powers disappear completely when 
P = 1, which is the case of interest for calculating the order parameters. 

The above results are presented in the hope that they may be a step towards 



verifying the conjecture (17) for the chiral Potts model order parameters. 
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